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ABSTRACT 

Recently, many efforts have been made to build dark energy models whose equation-of-state 
parameter can cross the so-called phantom divide Wde = — 1- One of them is the so-called hessence 
dark energy model in which the role of dark energy is played by a non-canonical complex scalar 
field. In this work, we develop a simple method based on Hubble parameter H (z) to reconstruct the 
hessence dark energy. As examples, we use two familiar parameterizations for H(z) and fit them to 
the latest 182 type la supernovae Gold dataset. In the reconstruction, measurement errors are fully 
considered. 
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I. INTRODUCTION 



Dark energy [1J has been one of the most active fields in modern cosmology since the discovery of 
accelerated expansion of our universe @, S, II i, S, 0, EH ■ The simplest candidate of dark energy is a tiny 
positive cosmological constant. However, as well-known, it is plagued with the "cosmological constant 
problem" and "coincidence problem" [l[ . In the observational cosmology of dark energy, equation-of-state 
parameter (EoS) Wde = Pde/ Pde plays an important role, where pd e and pde are the pressure and energy 
density of dark energy respectively. The most important difference between cosmological constant and 
dynamical scalar fields is that the EoS of the former is always a constant, —1, while the EoS of the latter 
can be variable during the evolution of the universe. 

Recently, evidence for Wd e (z) < — 1 at redshift z < 0.2 ~ 0.3 has been found by fitting observational 
data (see [g, d, E3, EH, EH EH, EH, EE EB] for examples). In addition, many best-fits of the present value of 
Wde are less than —1 in various data fittings with different parameterizations (see [I?} for a recent review). 
The present data seem to slightly favor an evolving dark energy with Wde being below —1 around present 
epoch from Wde > — 1 in the near past @. Obviously, the EoS cannot cross the so-called phantom divide 
Wde — — 1 for quintessence or phantom alone. Some efforts have been made to build dark en ergy model 
whose EoS can cross the phantom divide (see for examples 0, El, El 0, HH HI 0, HI [H, Ha H3, HI, 
[&5L I66L Wf\ and references therein) . 

In [9(, Feng, Wang and Zhang proposed a so-called quintom model which is a hybrid of quintessence and 
phantom (thus the name quintom). Phenomenologically, one may consider a Lagrangian density dHU, H3| 

Cquintom = ^ (d^lf ~ 7J O^s) 2 - V ((f>i , <f> 2 ) , (1) 

where <pi and 4>2 are two real scalar fields and play the roles of quintessence and phantom respectively. 
Considering a spatially flat Friedmann-Robertson- Walker (FRW) universe and assuming the scalar fields 
4>\ and 02 are homogeneous, one obtains the effective pressure and energy density for the quintom, i.e. 

Pquintom = ~ \^<t>\ ~ 7(01, 2 ), Pquintom = ^01 ~ + V {4>X, 0a), (2) 

respectively. The corresponding effective EoS is given by 

<Pj-4-2V(<f> 1 ^ 2 ) 
4^-^ + 2V(<k,(h) 



pquintom — - 9 - 9 „ xr / \ ' V / 



It is easy to see that w qu i n tom > — 1 when <p\ > §\ while w qu i n tom < — 1 when 4>\ < 4>\. The transition 
occurs when cf)\ = The cosmological evolution of the quintom dark energy was studied in [23|, l2~4j . 
Perturbations of the quintom dark energy were investigated in [29L [301 ] ; and it is found that the quintom 
model is stable when EoS crosses —1, in contrast to many dark energy models whose EoS can cross the 
phantom divide [28j ] . 

In E§| . by a new view of quintom dark energy, one of us (H.W.) and his collaborators proposed a 
novel non-canonical complex scalar field, which was named "hessence" , to play the role of quintom. In 
the hessence model, the phantom-like role is played by the so-called internal motion 8, where 9 is the 
internal degree of freedom of hessence. The transition from Wh > — 1 to Wh < — 1 or vice versa is also 
possible in the hessence model [18| . We will briefly review the main points of hessence model in Sec. [TTJ 
The cosmological evolution of the hessence dark energy was studied in [l9| and then was extended to the 
more general cases in [2(| ■ The w-w' analysis of hessence dark energy was performed in [2l| . 

In this work, we are interested in reconstructing the hessence dark energy. In fact, reconstruction of 
cosmological models is an important task of modern cosmology. For instance, the inflaton potential recon- 
struction was extensively studied in [3l| and references therein. The parameterizations and reconstruction 
of quintessence /phantom was considered in [H, HH . The other recent reconstructions of quintessence also 
include e.g. 34], HH, [H, H3, [H, [HH^I. The reconstruction of k-essence was studied in [4l|, HH, l4p. For 
the reconstructions of other cosmological models, see 0, |4|| [62|, [68|, H^] for examples. We refer to |46l| for 
a recent review on the reconstructions of dark energy. In this paper, after a brief review of the hessence 
model, we develop a simple method based on the Hubble parameter H(z) to reconstruct the hessence 
dark energy. As examples, we use two familiarparameterizations for H{z) and fit them to the latest 
182 type la supernovae (SNe la) Gold dataset [7J. In the reconstruction, measurement errors are fully 
considered. 
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II. HESSENCE DARK ENERGY 

Following [l8|, [l9j , we consider a non-canonical complex scalar field as the dark energy, namely hessence, 

$ = 01+^02, (4) 

with a Lagrangian density 

c h = \[ {d^f + {d^*f } - u(<p 2 + <r 2 ) = I [ {d^f - ^{d,ef } - v(<t>), (5) 

where we have introduced two new variables (0, 9) to describe the hessence, i.e. 

01 = 0cosh0, 2 = 0sinh6>, (6) 

which are defined by 

2 = 2 - 2 , coth0=^i. (7) 

02 

In fact, it is easy to see that the hessence can be regarded as a special case of quintom dark energy in 
terms of 0i and 02. Considering a spatially flat FRW universe with scale factor a(t) and assuming and 
9 are homogeneous, from Eq. we obtain the equations of motion for and 9, 

+ 3H 4> + 4>9 2 + y = o, (8) 

<f) 2 9 + (200 + 3H<t> 2 )9 = 0, (9) 

where H = a/a is the Hubble parameter, a dot and the subscript ", 0" denote the derivatives with respect 
to cosmic time t and 0, respectively. The pressure and energy density of the hessence are 

Ph = \(^- - y W). Ph = \{^~ + V(4>), ( 10 ) 

respectively. Eq. © implies 

Q = a 3 <p 2 9 = const. (11) 

which is associated with the total conserved charge within the physical volume due to the internal 
symmetry [18j, llfij . It turns out 

A Q 



a 3 2 

Substituting into Eqs. |SJ and (fTT)|) . they can be rewritten as 



(12) 



+ 3H0+^ + y = O, (13) 

to=\p-S?- y M> n=\p-Sp +y M- (i4) 

It is worth noting that Eq. (|13|l is equivalent to the energy conservation equation of hessence, namely, 
ph + SH (ph + Ph) = 0. The Friedmann equation and Raychaudhuri equation are given by 



" 2 = 3M^ iph + Pm) ' (15) 



H ' = —^772 (Ph + Pm + Ph) , (16) 



2M£ 



where p m is the energy density of dust matter; M p i = (8-kG)^ 1 ^ 2 is the reduced Planck mass. The EoS 
of hessence Wh = Ph/ Ph- It is easy to see that Wh > — 1 when 2 > <3 2 /(a 6 2 ), while Wh < — 1 when 
2 < Q 2 /(a 6 2 )- The transition occurs when 2 = Q 2 /(a 6 2 ). We refer to the original papers fl8l Il9j 
for more details. 



4 



III. RECONSTRUCTION OF HESSENCE DARK ENERGY 

Here, we develop a simple reconstruction method based on the Hubble parameter H (z) for hessence 
dark energy. From Eqs. (|15p and (fTH]) . we get 

V(cj>) = ZM 2 pl H 2 + M 2 pl H - l - Pm , (17) 

and 

= -2M 2 H-p m . (18) 



a 6 (f> 2 ~" pl 
Note that 

az 

for any function /, where z = a^ 1 — 1 is the redshift (we set ao=l; the subscript "0" indicates the present 
value of the corresponding quantity). We can recast Eqs. (fT7| and (fT5|) as 

dH 1 

= 3Af 2 tf 2 - M 2 (1 + z)ff— - -p m0 (l + z) 3 , (19) 

- |J(1 + z) 4 ^ 2 = 2M 2 pl (l + zY'H- 1 ^ - p m0 (l + z)H- 2 . (20) 
Introducing the following dimensionlcss quantities 

f>_ y 7 ^ ^ /S- g , 91 x 

F =^P!' ° = ^PV (21) 

Eqs. (fT9|) and ([20|) can be rewritten as 

V(z) = 3H 2 - (1 + z)H^- - ^ m0 (l + z) 3 , (22) 
dz 2 

(if) -Q 2 ^ 2 ( 1 + z ) 4 ^ 2 = 2 ( 1 + z ) _1 ^- 3 ^»o(l + ^- 2 , (23) 

where fl m o = Pmo/(3Af 2 ; iJQ ) is the present fractional energy density of dust matter. Once the H(z), or 
the H(z), is given, we can reconstruct V(z) and <j>(z) by using Eqs. (f2"2"|) and (|2"3"|) respectively. Then, the 
potential V{<j>) can be reconstructed from V{z) and 4>{z) readily. By using Eqs. (fTS|) and ([2"2"|l . we can 
reconstruct the EoS of hessence 

«*« - ^ - ' 1 + Ia + Z)L ^. (24) 

^ i-ano^- 2 (i + ^) 3 

The deceleration parameter 

*)--^=- i -i=- i+ < i+ ^- , f- < 25 > 

After all, it is of interest to reconstruct the kinetic energy term of hessence, K = 4> 2 /2 — Q 2 / (2a 6 (f) 2 ). 
From Eq. (fT8| . we have 

K(z) = = (1 + z)H^- - hl m0 {l + z) 3 . (26) 

It is worth noting that the reconstruction method presented here is sufficiently versatile for any H(z). 
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IV. EXAMPLES 



In this section, as examples, we consider two familiar parameterizations for H (z) and fit them to the 
latest 182 SNe la Gold dataset Q. And then, we reconstruct the EoS of hessence Wh{z), deceleration 
parameter q(z), the kinetic energy term of hessence K(z), the potential of hessence V(z), and the cj)(z) 
as functions of the redshift z. Also, we reconstruct the potential of hessence as function of cf>, namely 
V(<f>). In our reconstruction, measurement errors are fully considered. 



A. Parameterizations for H(z) and the latest 182 SNe la Gold dataset 

The latest 182 SNe la Gold dataset compiled in [7j provides the apparent magnitude m(z) of the super- 
novae at peak brightness after implementing corrections for galactic extinction, K-correction, and light 
curve width-luminosity correction. The resulting apparent magnitude m(z) is related to the luminosity 
distance (Il{z) through (see e.g. [5a |) 

m th (z) = M(M,H ) + Mog w D L (z), (27) 

where 

D ^ = ^^[ d ~ Z H { zJrLeters) < 28 > 

is the Hubble-free luminosity distance HqCLl/c in a spatially flat FRW universe (c is the speed of light); 
and 

M = M + 5 log 10 (^^) +25 = M — 5 log 10 h + 42.38 (29) 

is the magnitude zero offset (h is Hq in units of 100 km/s/Mpc); the absolute magnitude M is assumed 
to be constant after the corrections mentioned above. The data points of the latest 182 SNe la Gold 
dataset compiled in 0] are given in terms of the distance modulus 

Hobs(zi) = m obs (zi) - M. (30) 

On the other hand, the theoretical distance modulus is defined as 

fHh(zi) = m th (zi) - M = 51og 10 D L (zi) + fJ, , (31) 

where 

Mo = 42.38 - 51og 10 /i. (32) 
The theoretical model parameters are determined by minimizing 

X^parameters) = V (33) 
i=l ° ( Zl > 

where a is the corresponding la error. The parameter /zo is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over /io- However, there is an alternative 
way. Following (58l . l59l . l64j | , the minimization with respect to /io can be made by expanding the \ 2 °f 
Eq. (f3"3"| with respect to /io as 

\ 2 (parameters) = A — 2(1qB + (IqC, (34) 

where 

A(parameters) - V" [^obs(^) ~ r n th (zj] ^ = 0, parameters)] 2 



a m nh , ( z i) 



G 



^ m obs (zi) - m t h(zf,fi = 0, parameters) 
^(parameters) — >^ ~ 7 — n : 



LN2 



Eq. (|34|) has a minimum for /iq = -B/C at 



\ 2 (parameter s) — A(parameters) 



B (parameters)'' 
C 



(35) 



Therefore, we can instead minimize x 2 which is independent of /uq, since Xmin = Xmin obviously. 








-4 
Ai 







FIG. 1: The 68% and 95% confidence level contours in the A1-A2 parameter space for Ansatz I with the prior 
O to q = 0.30. The best fit parameters are also indicated by a solid point. 



In this work, we consider two familiar parameterizations for H(z) and fit them to the latest 182 SNe la 
Gold data 0]. At first, we consider the Ansatz I with 



H(z) = H [O m0 (l + zf + Ai(l + z) + A 2 (l + zf + (1 - Q m0 - A 1 - A 2 )] 



1/2 



(36) 



which has been discussed in [12, [H, \M, \M, [62, mm- Obviously it includes ACDM and XCDM with 
particular time- independent EoS of dark energy as special cases. As shown in [4l|, HtJ, even for the 
cases where this ansatz is not exact, one can recover the luminosity distance to within 0.5% accuracy 
using this ansatz in the relevant redshift range for the old 157 SNe la Gold dataset @- Therefore, this 
ansatz is trustworthy to some extent. Here, by fitting it to the latest 182 SNe la Gold data 0], for 
the prior tt m o = 0.30 [721 ]. we find that the best fit parameters (with la errors) are A\ — —3.28 ± 2.11 
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and A2 = 1.36 ± 0.84, while Xmm = 156.53 for 180 degrees of freedom. The corresponding covariance 
matrix [6(| (see also |47|) is given by 



Cov{A 1 ,A 2 ) = 



4.469 -1.774 
-1.774 0.711 



(37) 



In Fig. [TJ we present the corresponding 68% and 95% confidence level (c.l.) contours in the A\-A 2 
parameter space. 








-2.25 -2 -1.75 -1.5 -1.25 -1 -0.75 

w 

FIG. 2: The 68% and 95% c.l. contours in the wo-wi parameter space for Ansatz II with the prior Q m o = 0.30. 
The best fit parameters are also indicated by a solid point. 



Next, we consider the Ansatz II with 



H(z) = H a 



Q m0 (l + zf + (1 - fi m0 ) (1 + z)^ 1+Wo+w ^ exp 



3wiz 
l + z 



1/2 



(38) 



which is in fact equivalent to the familiar parameterization w c i e — wq + w\z/ '(1 + z) 3, EH, l6lL [H, r70l |. 
By fitting it to the latest 182 SNe la Gold dataset 0], for the prior Q m0 = 0.30 [72], we find that the 
best fit parameters (with la errors) are wo = —1.43 ± 0.32 and W\ — 2.79 ± 1.55, while Xmin = 156.56 
for 180 degrees of freedom. The corresponding covariance matrix reads 



, 0.101 -0.466 

Cov(w , Wl )=\ _ QAm 2AQ7 



(39) 



In Fig. [21 we present the corresponding 68% and 95% c.l. contours in the wq-w 1 parameter space. 
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FIG. 3: The reconstructed best fit Wh(z), q(z), V{z) and K (z) (thick solid lines) with the corresponding la errors 
(shaded region) for Ansatz I. 

B. Reconstruction results 

In our reconstruction, measurement errors are fully considered. The well-known error propagation 
equation for any y(x\, X2, ■ ■ ■ , x n ), 

^) = E£ cov( Xi ,x i )+2j2 £ (It It) (40) 

i \ 1 / x=x 1=13=1+1 \ 1 3 / x=x 

is used extensively (see 60] for instance). For Ansatz I with the prior il m o = 0.30, by using Eqs. (JUJ), 
(|24 |) -([26 |) . (|37|) . (|40|) and the corresponding best fit values of A\ and A2, we can reconstruct the EoS of 
hessence Wh(z), deceleration parameter q(z), the potential of hessence V(z), and the kinetic energy term 
of hessence K(z) as functions of the redshift z, with the corresponding ler errors. We show the results 
in Fig. O It is easy to see that crossed —1 and the universe transited from deceleration (q > 0) 
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to acceleration (q < 0); the reconstructed Wh(z) is well consistent with the three uncorrected Wo. 25, 
Wo. 70 and W1.35 data points with their corresponding ler error bars for the weak prior Q which are 
model-independent. 

However, the error propagation equation (I40p is invalid when we reconstruct the <fi( z ) an d hence the 
V ((/)), since 4>{z) is obtained from a differential equation, i.e. Eq. P3"j). To evaluate the error propagations, 
we use the Monte Carlo method instead. That is, we generate a multivariate Gaussian distribution from 
the best fit parameters and the corresponding covariance matrix. And then, we randomly sample N 
pairs of the parameters {^1,^2} from this distribution. For each pair of {Ai,^}, we can find the 
corresponding <p(z) and V(z) from Eqs. (f/!3"|) and (J22) respectively. Hence, the V{4>) is in hand. Finally, 
we can determine the mean and the corresponding la error for the 4>{z) and V(<p) from these N samples, 
respectively. In Fig. HI we show the reconstructed (j>(z) and V((f>) with the corresponding la errors for 
Ansatz I with the prior Q m Q = 0.30. In which, we have used the demonstrative initial value <po = 0.05 at 
z = and Q = 1, and have chosen the solution with d<j)/dz > for the reconstructed (j>(z); we have done 
N = 1000 samplings. 



O 
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FIG. 4: The reconstructed (/>(z) and V(<j>) (thick solid lines) with the corresponding la errors (shaded region) for 
Ansatz I. See text for details. 



For Ansatz II, the method to reconstruct the EoS of hessence Wh(z), deceleration parameter q(z), the 
kinetic energy term of hessence K (z), the potential of hessence V(z), the </>(z) as functions of the redshift 
z, and the potential of hessence as function of <fr, namely V(<p), is the same for Ansatz I. We present the 
results in Figs. [5] and [U Once again, it is easy to see that Wh crossed —1 and the universe transited from 
deceleration (q > 0) to acceleration {q < 0); the reconstructed Wh{z) is well consistent with the three 
uncorrelated Wo. 25, Wo. 70 and W1.35 data points with their corresponding la error bars for the weak 
prior Q which are model- independent. 
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FIG. 5: The reconstructed best fit Wh(z), q(z), V(z) and K(z) (thick solid lines) with the corresponding ler errors 
(shaded region) for Ansatz II. 

V. DISCUSSION 

In this work, we have developed a simple method based on the Hubble parameter H (z) to reconstruct 
the hessence dark energy. If the observational H(z) is obtained, the reconstruction of hessence dark 
energy is ready. It is worth noting that the reconstruction method presented here is sufficiently versatile 
for any H{z). 

As examples, we reconstructed the hessence dark energy with two familiar parameterizations for H(z). 
It is easy to see that this reconstruction method works well. However, these parameterizations for H (z) 
are not the direct measure of H (z) from observational data. We can just say that they are consistent with 
the recent observational data. Can we be in a more comfortable situation? In fact, some efforts are aiming 
to a direct measure of H(z) from observational data. Analogous to the estimates of w(z) [8], a method 
to obtaining the estimates of H(z) is proposed in [ll| (see also Actually, in the latest paper Q by 

the Supernova Search Team led by Riess, a rough estimate of H(z) is obtained by using the new Hubble 
Space Telescope discoveries of SNe la at z > 1. Other new method to determine the Hubble parameter 



11 



as function of redshift, H(z), is also proposed in [5(| recently. In a very different way, by using the 
differential ages of passively evolving galaxies determined from the Gemini Deep Deep Survey (GDDS) 1511 
and archival data [52j |. Simon et al. determined H(z) in the range o<z<i.8 0, 0, HM lij. 
However, up to now, all observational H(z) obtained by various methods are too rough to give a reliable 
reconstruction of hessence dark energy. A good news from [55] is that a large amount of H(z) data 
is expected to become available in the next few years. These include data from the AGN and Galaxy 
Survey (AGES) and the Atacama Cosmology Telescope (ACT), and by 2009 an order of magnitude 
increase in H(z) data is anticipated. Therefore, we are optimistic to the feasibility of the reconstruction 
method of hessence dark energy proposed in this work. 
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FIG. 6: The reconstructed (j)(z) and V{<j>) (thick solid lines) with the corresponding la errors (shaded region) for 
Ansatz II. See text for details. 
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